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ABSTRACT 

As a neutron star spins down, its core density increase, changing the relative equi- 
hbrium concentration, and causing deconfinement phase transition as well. Hadron 
matter are converted into quark matter in the interior, which enhances the deviation 
of chemical equilibrium state. We study such deviations and its chemical energy re- 
lease. Applying to the simulation of cooling neutron stars, we find the surface effective 
temperature of neutron stars is promoted obviously. This implies that the deconfine- 
ment phase transition is able to raise the chemical heating efficiency. 
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1 INTRODUCTION 

Neutron star cooling is an important tool to probe 
its inner structure. The comparison of cooling models 
with observation of thermal emission provide insight into 
the EO S of dense matte r, and the signature s of e xotic 
particles (|Page et al.l bood l (|Yakovlev fc PethicS |2004| ). A 
newly formed neutron star loses the thermal energy through 
neutrino emission initially, and is taken over by the surface 
photon radiation about ~ lO'^yr after birth. However, heat- 
ing sources inside neutron stars can keep neutron stars hot 
beyond the standard cooling timescale ~ f O^yr. Some mech- 
anisms have been extensively discussed. These include the 
mutual fric tion between superfluid and normal component 
of the star (|Shibazaki fc Lamblll989l). heat released t hrough 
deconfinement phase transition! Kang fc Zhen3l2007l )and re- 
lease of strain energ y stored b y the solid crust due to spin- 
down deformation (j Cheng et al. 1992). 

Another of these mechanis m is rotochemical heating, 
firstly proposed by Reiseneggerl Reiseneggei 19951) and ap- 
plied to quark stars by Cheng! Cheng fc Pail " 199^ ). Then 
Reisenegger improved it, in the framework of general rel- 
ativity, b5'_consid^ring_tlM_J^ via realistic 
EOSs JPernandez fc ReisenggeilbOOSl) and superfluid nucle- 
ons core( Petrovich fc ReisenggedbOld ). All studies above 
agree that rotochemical heating is important for old neu- 
tron stars. 

Rotochemical heating origins in the deviation from 
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weak reaction equilibrium due to spin-down compression. 
As the star spins down, the decrease of the centrifugal force 
leads to compression, increasing the star's internal density. 
This compression results in a displacement of the equilib- 
rium concentration of each particle species and changes the 
chemical equilibrium state. Reactions that change the chem- 
ical composition drive the system to the new equilibrium 
configuration. If the reaction rate is slower than the change 
of the equilibrium concentrations due to spin-down com- 
pression, the star is never exactly in chemical equilibrium, 
so energy is stored. The excess of energy is dissipated by 
enhanced neutrino emission and heat generation. 

It has long been hypothesized that some compact stars 
might actually be "hybrid stars" containing cores of quark 
matter (|Bavm fc Chin|[l973 ). Unlike neutron stars with pure 
hadronic matter and strange stars with pure strange mat- 
ter, those hybrid stars probably contain the m ixed phase in 
light of Glendenning's phase transition theorv l|GIendenning| 
I2000I ). Of course, the compression of a hybrid star due to 
spin-down causes deconfinement phase transit ion (hereafter 
DPT), besides nonequilibium Urea processes tSteiner et al.l 
2009). The spin-down compression increases the interior 
density gradually, which results in the transformation of 
hadron matter into quark matter in the interior little by lit- 
tle. The deconfinement reactions change relative concentra- 
tions of the particles, neutrons, protons, and u, d, s quarks, 
and result in the displacement of equilibrium state as well. 
Because neutrons and protons are converted into two fla- 
vor quarks, and it takes a finite amount of time that two 
flavor quarks decay into three. Thus the deconfinement con- 



2 Wei Wei Xiao-Ping Zheng 



state 1 state2 

KN^ \aaC 

K N, N. 




Figure 1. The DPT lead to the change of quark and hadron 
numbers, which was treated in equilibrium condition from equi- 
librium state 1 to equilibrium state 3 in previous study, AA^q = 
AAf^ + AN'^ + AN!,, ANh = AN^ + AN^. In fact, neutrons and 
protons are converted into u and d quarks, without s quarks. 
So the deconfinement reactions break the equilibrium state and 
drive the system deviating to a nonequilibrium state 2, ANq = 
ANu + AN^,ANh = AN„ + ANp. Weak reactions adjust the 
concentrations of these particles and try to recover the system to 
a new equilibrium state 3. 



tribution need to be considered in the calculation of the 
/3-equilibrium state deviation. 

The purpose of this paper is to study the effect of DPT 
on rotochemical heating in mixed phase. In section 2, we at- 
tempt a model to describe the deviations from /3-equilibrium 
state under the circumstance of DPT. In section 3, we de- 
rive the equations of chemical and thermal evolution, and 
discuss an application of a uniform star model with mixed 
phase matter. The results and corresponding explanations 
are presented in section 4. In section 5, we give our conclu- 
sions and suggest future improvement. 



2 MODEL 

Hybrid stars have more complicated matter composition 
than purely neutron stars, which contains neutrons, protons, 
u, d and s quarks and electrons. The relative concentrations 
of these particles are adjusted by the following weak reac- 
tions. 



It -I- e~ d + Ve, 
u -\- e~ s + Ue. 
p + e~ n + Ue 



d — !> It + e 

s u + + Ve, 

n ^ p + e~ + Ve. 



(1) 



At the same time, the compression leads to the deconfine- 
ment reactions. 



2d + u, p ^ 2u + d. 



(2) 



Nucleons dissolve into only two flavor quarks, so the relative 
concentrations of three flavor quarks need to be adjusted by 
non-leptonic reaction. 



ti + d o s - 



(3) 



Thus there are two channels, through which the /3- 
equilibrium is violated. One is a direct contrib ution of spin- 
down compression discovered by Reisenegger (jReiseneggerl 



Il995ll . Another is the contribution of deconfinement reac- 
tions caused by the spin-down. The deconfinement reactions 
take place firstly and then result in a departure from /3- 
equilibrium state. 

The deconfinement reactions lead to the change of quark 
baryon number density Auq and change of nucleon baryon 
number density Auh- Since neutrons and protons are con- 
verted into u and d quarks, without s quarks, we have 



An, = Auu + Arid, 



Auh = Ann + Aup. 



(4) 



The reactions© change the concentration of particles and 
drive the system deviating from an equilibrium state (statel) 
to a nonequilibrium state (state 2). Then, the system tries 
to recover to a new equilibrium state (stateS) through re- 
actions Qand ((31 . We give a description of these processes 
in Figure 1. Since the deconfinement reactions is faster than 
weak reactions, the system is never exactly in chemical equi- 
librium due to deconfinement behavior. 

For quarks, the non-leptonic reaction Q and Durca 
reactions ([1} are all involved in the evolution of chemical 
potential, among which only two processes are needed in 
the calculation of the chemical departure, resulting in roto- 
chemical heating. In the later calculation, we find that the 
contribution of the nucleon dominates hybrid star thermal 
evolution and the chemical evolution of nucleon is indepen- 
dent of reaction So we don't distinguish the difference 
between the non-leptonic reaction and Durca reactions, and 
choose the Durca reactions to calculate the chemical poten- 
tial difference. 

To express the chemical deviation quantities, under con- 
sideration of deconfinement reactions, we define the chemi- 
cal potential of particle i in nonequilibrium state (state2)as 
fi"'' and Hi in new equilibrium state (state 3) . The chemical 
potential differences of three channels in reactions(IT| are: 



Shd 
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- Sfln 


- Sue, 


SfJ,s 


= Sfls 
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SfiN 


= S/In 


- S^p 
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where reaction channels denoted by capitals, particles de- 
noted by small letters. 

We assume n„ = anh,np = {l~a)nh , under beta equi- 
librium, then An,! — aAuh, Aup = (1 — a)Anh, this leads 
to the following equality, with constraint by reactions ([2}, 



"+ 1 A 

And = t; An„, 

2 — a 

Substitute ([5]) into (Q, we obtain 
2 - a 



(5) 
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At zero temperature, chemical potential of Fermions satisfies 
j-ij — voj + (37r^n/) 3 , then we determine chemical potential 
of particles in each state for density change and density of 
a given initial state, 
state 1: 

^.u = (37r ng)3(3+ — -) 3, 

A*d = Ms = Mil + Mei 

, 2 a.„^3/ie 3m^ -f3m^ 1 
Atp = (37r nh) 3(2 -f — ) 3, 
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state 3: 





2 — a An„ 



1 + a An, , 
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Mn - ■m'i An„ 



Mu = f^u[l + - — ), 

I /-I I '^^9 \ /I I A'^9 \ I ' ' I ' 

Md = Md(l + ^5 ) - Me(l + ^ ) + Me = Mu + Me, 

6nq 6nq 

I / 

Ms — Md, 

Mn Mp Me ■ 

State 2 and 3 can be expressed in light of state 1 and density 



change, we immediately write 



Sjiu = m" 



Mu 



f no / c- no 



nq 
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c nq I c nq I r nq I 

0[l,n= iln ~ [i-n-, O/lp = flp^ ~ flp, d^e^f^e " Me • 

since state 3 has been defined as equilibrium state, the /3- 
equilibrium condition should enter to simplify our calcula- 
tions. So 
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3 EQUATIONS OF CHEMICAL AND 
THERMAL EVOLUTIONS 

According to the discussion above. We first calculate the 
departure rate from /3 equilibrium in channel d u+e^ -\-Vc 

dSflJ^ 1 _ ^q y^q nq-\ 
; = lim Mj — M7/ Me 

To leading order, we obtain 
^dS^D 1 _ 2a — 1 /id duq ^ ^ 



dt 



\departure 



3 Uq dt 



In the similar way, we give 



.d^Ms 
^ dt 

,dSfAN_ 
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,0 2 a-, 1 rfrih 



({3n^np)3 + 

On the other hand, the departure is preferred to recover its 
equilibrium through reactions U]). Therefore, we have 

rdSnd. 



dt 



[ ]restore — Ts (/i) , 



^ dt 



= -r„E"^{fi), 



where x indicates d. 
sity, Exx is partial 



s, n fraction of baryon den- 

de rivati ve of chemical energy 

per barvon (|Cheng fc Pail 1 19961 ). In the hadron-quark 
mixed p hase, hadron and q uark undergo direct Urea 
processes teang fc Zhene 2007 ). Fd, F aand r„ are th e reac- 



tion rates per barvon (|Dai et al.lll993l ) (|Haensellll993 ). Com- 
bine departure and corresponding restore, we immediately 
get the time evolution equations of the chemical potential 
differences. 

dSfiD 



dt 



dSfis 
dt 



dSfj,]\ 
dt 



^dExxil^), 



■ a M< 



dn„ 



dt 



[Siv^np) 



((37r2n„)3 +1711)2 ((37r2np)3 -l-m^); 



(6) 



(7) 



(8) 



Equation ((6])and (O show that D and S channels obey their 
evol utional equations re spectively, differing from Cheng & 
Dai (|Cheng fc Daill 19961 ) discussion of strange quark matter. 
Inclusion of deconfinement effect causes such serious differ- 
ence. To solve these equations, thermal equation is neces- 
sary, reads 

dT 

Cv—r- = (Fd^Mo+rs^MS + F„(5/ijv) , ^ 

dt (9) 

- (ed -I- ts + En) — -B7, 

where the three terms represent the chemical energy re- 
leased by the reactions, the energy radiated by neutrino and 
antineutrino, and the energy in photons from the surface 
of the star. c„ is the specific heat of dense matter c„ — 
XCq + (1 — x)ch- Now we need to treat fj,f,nq,nh as the func- 
tion of total baryon number density. In 1992, Glendenning 
first presented the existence of first -order phase transition 
with Gibbs construction. He applied the relativistic mean 
field theory for hadronic matter and MIT for quark matter to 
construct the phase diagram. Baryon densities of each phase 
in m ixed phase are depicted in his figure? (|Glendenning| 
Il992h . We here consider a uniform star model for simula- 
tion and take baryon density n — 1.0 fm~^ .We approxi- 
mate, around the baryon density, n, « 0.3431 -I- 0.8228n, 
rih ^ 0.07143 + 0. 7428n, Md ~ jMn « |(4.826 + 1.44n). 
From Reisenegger(Rcise neggeill 19951 ). the total baryon num- 
ber density changes with the spin-down by 

1 dn Q dfi 

n dt Gpc dt 

where pe is central density of the stars. 



4 RESULTS 



1ms and a = |, the 



We give an initial rotation period Pi 
latter corresp onding to the thresho ld for direct Urea process 
to be allowed iLattimer et al.|[l991^ . After giving the strong 
coupling constant, s-quark mass, and the surface magnetic 
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Figure 2. Chemical evolution. The variables S^/-Kk{a, measure for 
the dcpature from chemical equilibrium) is plotted logarithmically 
as functions of time. The solid line, dashed line and dotted line are 
chemical evolution for S channel, D channel and nucleon channel 
respectively. The four panels correspond to different magnetic field 
strengths: (a)B = IQ^G; (b) B = lO^G; {c)B = lO^^G; (d) B = 
10"G. 



field strength , we integrate tlie coupling evolution equations 

Figure 2 sliows tiie ciiemical evolution in liadron-quark 
mixed piiase. The solid line, dashed line and dotted line 
are chemical evolution for S channel, D channel and nu- 
cleon channel, respectively. Since neutrons and protons are 
converted into u and d quarks, without s quarks, the chem- 
ical potential departures of S channel is higher than that 
of D channel. The four panels correspond to four differ- 
ent magnetic field strengths:(a)B = l(fG, (b) B = lO^G, 
{c)B = 10^°G, (d) B = 10"G. The stronger the magnetic 
field strength, the more the kinetic energy of rotation con- 
verting into chemical energy ,the larger the chemical poten- 
tial departure, the earlier the appearance of heating effect. 
The platform indicates a relative quasi-equilibrium process. 
And the rapid dropping of chemical departure weaken the 
rotochemical heating at later times for strong magnetic. But 
rotochemical heating is still substantial for B — 10* G at 
later times. Figure 3 gives the chemical potential departure 
of nucleon channel with magnetic field strengths B = 10*G 
and B = 10^° G, to compare deconfinement and spin-down 
situations. Obviously, the deconfinement reactions enlarge 
the chemical potential departure. 

The chemical evolution infiuence the cooling of neutron 
stars, when the chemical heating mechanism is considered. 
The DPT enhance the chemical potential departure of par- 
ticles, and effect the thermal evolution of stars containing 
mixed phase matter as well. Figure 4 shows the evolution 
of the effective surface temperature for uniform hybrid star 
with rotochemical heating under consideration of DPT, and 
the traditional case. Although we use a toy model in which 
several approximations are made, the effective surface tem- 
perature of the star is about 40 percent higher than the 
traditional one, and the heating effect appears earlier. 

We also analysis the contributions of quarks and nucle- 
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Figure 3. Chemical potential departure is plotted logarithmically 
as a function of time for considering DPT {solid line )and without 
it{dashed line ), with magnetic field strengths B = 10*G and 
B = IQi^G. 
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Figure 4. Effective surface temperature as a function of time, for 
rotochemical heating considering DPT {solid line )and traditional 
rotochemical heating {dashed line ),with magnetic field strengths 
B = lO^G and B = 10^°G. 



ons by comparing temperature magnitudes. Figure 5 gives 
the results with B = 10*G. Since nucleon's curve is close to 
full one, the contribution of the nucleon dominates hybrid 
star thermal evolution. 

Figure 6 shows the evolution of the effective surface 
temperature with different magnetic field strengths B = 
l(fG{dashed line), 10'-^G{dotted line), 10^° G{dashed-dotted 
line), 10^^ G {dashed- dotted- dotted line), and in the absence 
of chemical heating. One sees that the effect of rotochemical 
heating with DPT can be substantial at later times for stars 
with weak fields and still important at earlier times for stars 
with strong fields. 
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Figure 5. Eflfective surface temperature as a function of time, 
with magnetic field strengths B = 10® G, for only considering 
quark's chemical deviation and emission(£ias/ie(i line ),only con- 
sidering nucleon's(doMe(i line ) and considering hoth{solid line 
)■ 




Figure 6. Effective surface temperature as a func- 
tion of time, with no heating (solid lines) or roto- 
chcmical heating with magnetic field strengths B = 
10^G{dashed line), 10^G{doUed line), 10^^ G {dashed- dotted line), 
and 10^^G{dashed-dotted-dotted line). 



5 CONCLUSIONS 

Besides the spin-down, an adjoint of DPT indeed causes the 
departure from /3-equilibriurn in mixed phase matter. This 
mechanism is investigated in this paper. The new mechanism 
makes a substantial promotion in effective surface temper- 
ature. We usually believe hybrid stars cool down rapidly 
following standard cooling theory. Our results show that a 
hybrid star should has much higher temperature when the 
chemical heating with DPT is considered. Although we con- 
sider a uniform density model to facilitate our understanding 
of the mechanism, the results are very close to the hybrid 
stars containing large mixed phase core if their mean density 
equals the assumed uniform density. 

In our calculations, we set a as a constant and only 
consider the contribution of the DPT in chemical deviation. 



The new mechanism makes a greater contribution than the 
spin-down one. So the direct contribution of spin-down com- 
pression is not important for rotochemical deviations in stars 
containing mixed phase matter and can be ignored in calcu- 
lations. 

This paper focus on studying a new mechanism in heat- 
ing generation inside hybrid stars, our consideration of ini- 
tial period Pi = 1ms is reasonable for millisecond pul- 
sars, but is too short for classical pulsars with higher mag- 
netic field strengths. So the more reasonable initial period 
is needed in the accurate calculations comparing with the 
observations. 

Accurate calculations arc worth being expected. Our 
future works would study the thermal evolution of hybrid 
star including DPT, taking the structure of star into 
account in the frame of general relativity, using realistic 
EOSs of dense matter. The results thus would compare 
with the observations. 
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